Weighted Quadrature in Krylov Methods
نویسنده
چکیده
The Krylov subspace approximation techniques described by Gallopoulos and Saad 2] for the numerical solution of parabolic partial diierential equations are extended. By combining the weighted quadrature methods of Lawson and Swayne 6] with Krylov subspace approximations, three major improvements are made. First, problems with time-dependent sources or boundary conditions may be solved more eeciently. Second, methods are derived which have the stability properties (such as A{stability) of the underlying rational approximation to the exponential function. Third, it is possible to present methods which are robust under space discretization reenement. In particular, a xed precision is essentially maintained for the same time integration method and for constant values of the parameters, when the spatial resolution is increased.
منابع مشابه
Spectral Methods for Time-dependent Variable-coefficient PDE Based on Block Gaussian Quadrature
Block Krylov subspace spectral (KSS) methods are a “best-of-both-worlds” compromise between explicit and implicit time-stepping methods for variable-coefficient PDE, in that they combine the efficiency of explicit methods and the stability of implicit methods, while also achieving spectral accuracy in space and high-order accuracy in time. Block KSS methods compute each Fourier coefficient of t...
متن کاملMatrices, Moments and Quadrature: Applications to Time-Dependent Partial Differential Equations
The numerical solution of a time-dependent PDE generally involves the solution of a stiff system of ODEs arising from spatial discretization of the PDE. There are many methods in the literature for solving such systems, such as exponential propagation iterative (EPI) methods, that rely on Krylov projection to compute matrix function-vector products. Unfortunately, as spatial resolution increase...
متن کاملEnhancement of Krylov Subspace Spectral Methods by Block Lanczos Iteration
Abstract. This paper presents a modification of Krylov subspace spectral (KSS) methods, which build on the work of Golub, Meurant and others, pertaining to moments and Gaussian quadrature to produce high-order accurate approximate solutions to variable-coefficient time-dependent PDEs. Whereas KSS methods currently use Lanczos iteration to compute the needed quadrature rules, our modification us...
متن کاملDerivation of High-order Spectral Methods for Time-dependent Pde Using Modified Moments∗
Abstract. This paper presents a reformulation of Krylov Subspace Spectral (KSS) Methods, which build on Gene Golub’s many contributions pertaining to moments and Gaussian quadrature, to produce high-order accurate approximate solutions to variable-coefficient time-dependent PDE. This reformulation serves two useful purposes. First, it more clearly illustrates the distinction between KSS methods...
متن کاملRecent Advances in Krylov Subspace Spectral Methods
This paper reviews the main properties, and most recent developments, of Krylov subspace spectral (KSS) methods for time-dependent variable-coefficient PDE. These methods use techniques developed by Golub and Meurant for approximating elements of functions of matrices by Gaussian quadrature in the spectral domain in order to achieve high-order accuracy in time and stability characteristic of im...
متن کامل